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Abstract: This paper explores various distributional aspects of random 
variables denned as the ratio of two independent positive random variables 
where one variable has an a stable law, for < a < 1, and the other 
variable has the law defined by polynomially tilting the density of an a 
stable random variable by a factor 8 > —a. When 8 = 0, these variables 
equate with the ratio investigated by Lamperti (66) which remarkably was 
shown to have a simple density. This variable arises in a variety of areas and 
gains importance from a close connection to the stable laws. This rationale, 
and connection to the PD(a, 8) distribution, motivates the investigations 
of its generalizations which we refer to as Lamperti type laws. We identify 
and exploit links to random variables that commonly appear in a variety 
of applications. Namely, Linnik, generalized Pareto and z— distributions. In 
each case we obtain new results that are of potential interest. As some high- 
lights, we then use these results to (i) obtain integral representations and 
other identities for a class of generalized Mittag-Leffler functions, (ii) iden- 
tify explicitly the Levy density of the semigroup of stable continuous state 
branching processes (CSBP) and hence corresponding limiting distributions 
derived in Slack, Zolotarev(87; 97l), which is related to the recent work by 
Berestycki, Berestycki and Schweinsberg , and Bertoin and LeGallQ 0) on 
beta coalescents. (iii) We obtain explicit results for the occupation time 
of generalized Bessel bridges and some interesting stochastic equations for 
PD(a, 0)-bridges. In particular we obtain the best known results for the 
density of the time spent positive of a Bessel bridge of dimension 2 — 2a. 
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1. Introduction 

Let S a , for < a < 1 denote a positive stable random variable having Laplace 
transform 




•Supported in part by grants HIA05/06.BM03, RGC-HKUST 6159/02P, DAG04/05.BM56 
and RGC-HKUST 600907 of the HKSAR. 



1 

imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. J ames / Lamperti Laws 



2 



Additionally, for 9 > —a define variables S a ,e independent of S a whose laws 
follow a polynomially tilted stable distribution having density proportional to 

t~ e f a {t)- When 9 = 0, S a . := S' a = S a . In this case Lamperti ((§§), (see also 



Zolotarev (|9q). and Chaumont and Yor (|2lf )). showed that, despite the general 
intractability of f a , the ratio 

v Si a 

has a remarkably simple density given as 

fx {y) = Sln(TO) - - ^ for „ > 0. (1.1) 

J " yy> 7T y 2a + 2y a cos(Tra) + 1 y v ' 

This variable arises in many important and often seemingly unrelated con- 
texts. For instance, (17; H; 27; 71;[z3 76|). Inspired by these facts and connec- 
tions to the (a, 6) Poisson Dirichlet family of distributions discussed in Pitman 
and Yor (|82l ) leads us to investigate properties of variables defined as 

v d S a 



We refer to these variables as being Lamperti variables or variables having Lam- 
perti type laws. Our purpose, from a broad perspective, is to demonstrate that 
these variables have strong connections to more familiar random variables that 
appear in a variety of applications in probability, statistics and related fields. 
In other words, the Lamperti variables, albeit often hidden, appear in many 
important contexts. Furthermore, we show how to utilize these links to both 
deduce properties of X a> $, and develop new non-trivial results related to the 
linked variables. These results can also be potentially used to expand modeling 
capabilities. Our results are suggestive of an active beta-gamma-stable calculus 
that extends the notion often associated with beta and gamma variables via 
Lukacsl69T) characterization. 



1.1. Outline 

We now present an outline of this paper highlighting specifics. More detailed 
references can be found in each section. We note that each section contains new 
results of a non-trivial nature that in some cases are generalizations of existing 
results. In addition, combined, they represent a nice partial survey of linked 
variables. Section 2 consists of essentially two parts. The first develops a series 
of pertinent distributional results for X a .e and for a more broader class de- 
fined by multiplying the Lamperti variables by beta variables. One shall notice 
the class of random variables we denote as X a a \ plays a major role through- 



out the sections. This multiplication is based on ideas we developed in (|45l ). 
The second constitutes a natural progression of ideas, each section building on 
the previous one. Specifically, section 2.3 establishes links with positive Linnik 
variables. Where, as highlights, we obtain expressions for the density of Linnik 
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variables and also establish an interesting gamma identity. Section 2.4, exploits 
this identity in connection with generalized Pareto distributions. Albeit brief, 
the main result is used to identify an unknown limiting distribution obtained by 
Zolotarev (f97h and Slack §8% which we discuss in section 4. Section 2.4, influ- 
enced by (| 14t 42; 84; 9^), uses the characterization in the previous sections to 
demonstrate how one can develop a calculus of sorts involving z-distributions. 
In particular, we identify new classes of random variables, arising as solutions to 
stochastic equations involving z-distributions, having both hyperbolic character- 
istic functions and variables whose density can be computed explicitly. Section 
3 obtains results for a generalization of Mittag-Leffler functions that can be ex- 
pressed as Laplace transforms of or X a ,g and can be represented in terms of 
densities of Linnik variables. Section 4 solves a fairly hard problem, identifying 
the explicit Levy density of the semigroup of stable continuous state branching 
processes. This section continues with the idea of using results from previous 
sections. Results in sections 5 and 6, with the exception of a = 1/2, present the 
best known results for occupation times of various quantities including times 
spent positive on certain random subsets. We also develop a series of interest- 
ing stochastic equations. As one highlight we obtain results for the otherwise 
elusive case of Bessel bridges of dimension 2 — 2a. For some related works see 
(HI 51; 52; 59|;[94|). Section 7 discusses aspects of Brownian time changed models 
where we close by exploiting an interesting, yet not well known, representation 
of symmetric stable variables of index < 2a < 1, found in (jjjl ). 



1.2. Some notation and background 

Here we briefly recount some notation and background related to Bessel pro- 
cesses and the Poisson Dirichlet family of laws. See Pitman (79|; 8^) for a more 



precise exposition. Let B :— (B t ,t > 0) denote a strong Markov process on M 
whose normalized ranked lengths of excursions, (Pi) E V = {s = (s%,S2,---) '■ 
si > S2 > ■ ■ ■ > and s « = ty> f° now a Poisson Dirichlet law with pa- 

rameters (a,0) for < a < 1, as discussed in Pitman and Yor Denote 
this law as PD(a, 0). We hereafter consider B := (B t , < t < 1). Let (L t ; t > 0) 
denote its local time starting at 0, and let Te = inf{£ : L t >£},£> denote its 
inverse local time. In this case r is an a-stable subordinator where we choose 
ti = S a . Writing L t = L(t), there is the scaling identity (see (|8lh ). 

L(t) s d_ Q 

t b~(s)j 



a ' 



where the inverse local time up to time 1, L±, satisfies 

Li := r(l - a)" 1 lim e a \{i : Pi > e}\ a.s. (1.2) 

e — >n 



and is said to follow a Mittag-Leffler distribution. This shows that (L t ,T t ) have 
distributions determined by PD(a,0). Furthermore, independent of (Pi), we 
suppose that for a fixed < p < 1, B is symmetrized so that P(B t > 0) = p. 
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Under these specifications B could be a p-skewed Bessel process of dimension 
2 — 2a. In particular when p = 1/2, a =1/2 then B behaves like a brownian 
motion. An interesting aspect of B is the time its spends on certain subsets of R. 
Let = J* * I( Bs>0 )ds and A[~^ — J * l( Bs<0 )ds, such that t = A^+A[~\ denote 
the time B spends positive and negative respectively up till time t. Remarkably, 
by excursion theory, the time changed-processes (A+;£ > 0) and (A^;£ > 0) 

are independent a-stable subordinators such that At = p l ^ a S a and A^ ^ = 

(1 - p) 1/a S' a , for S' a = S a . This leads to, 

v d d S a cX a a A (j , 

A Q = -^y = c— and — — = — = AJ. 1.3) 
AV, } S' a cX a + 1 t? 

Hereafter, denote the law that governs B and related functionals under the 
above specifications as P^f . Denote the corresponding expectation operator as 
E^q Thus writing V^^Af £ dx)/dx equates with the density of the time spent 



positive on [0, 1] of a p-skcwed Bessel process. As noticed by, ||3j ISJJ) , this law 



was originally obtained by Lamperti (|66f ) and from (|1.3p it equates with, 

P$>(^i e dx)/dx = P(cX a /(cX a + 1) e dsj/da;. 

Now for 6* > let P^g and E^g, denote the law and expectation operator 
of functionals connected to a p-skewed process whose excursion lengths follow 
PT)(a,9). In particular if (Pi) is distributed according to PD(a,9), then it sat- 
isfies, for measurable H, 

^m(Pi))} = ^^ E^[H(( Pi ))rr e ]. 
When 8 = a, this corresponds to the case of a Bessel bridge of dimension 2 — 2a. 

(br) 

We use the notation A\ for the variable that satisfies 

F^ e (A? r ^dx)=F^ e+a (A+edx). 

Note also that under P^fl, L\ = S~^, which is also equivalent in distribution to 
the a-diversity of a PD(a, 9) law. Let now Ui, U2, ■ ■ . , denote a sequence of iid 
Uniform[0, 1] random variables, for (Pi) distributed according to PD(a,9) and 
< u < 1, the class of PD(a,8) random cumulative distribution functions are 
defined as 

00 

P a ,e(u) = ^2P k I(u k <u)- 



Furthermore under P„ \. 



k=l 



A+ = P a , e (u) 



for a fixed u. See Bertoin (8j) for applications to coagulation/fragmentation 
phenomena where it is called a PD(a, 6*)-bridge and Ishwaran and James (|44r ) 
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[see also Pitman d23)]for applications to Bayesian statistics, where in particular 
P a fi is referred to as a Pitman- Yor process. Under this name the process has also 
been applied to problems arising in natural language processing, see Teh(89). 
When 9 > and a — Po,e is a Dirichlet process made popular by Ferguson (|34l ). 

For basic notation, we write j a and f3 a ^ to denote a gamma random variable 
with shape a and scale 1, and a beta random variable with parameters (a, b). 
Additionally £ CT denotes a Bernoulli variable with success parameter < a < 1 . 
If X and Y are random variables we will assume that XY is a product of 
independent random variables unless otherwise specified, if we write X, X' this 
will mean X = X' but they are not equal. Lastly we always consider c a = p/q = 
p/(l — p) where q = 1 — p unless otherwise specified. 



2. Distributional Results for X a ^g 



In this section we shall derive various distributional properties of X a ^g. For 
r > and < a < 1, we will sometimes work with the parametrization tct, to 
accommodate values such as tct = 9 > and tct = 9 + a. First, we briefly discuss 
some pertinent properties of random variables referred to as Dirichlet means and 
the related class of infinitely divisible random variables whose distributions are 
generalized gamma convolutions (GGC), as they will play a significant role in 
our exposition. For more details and related notions, one may consult earn 
2^h H3; [H; [III; 0) and, in particular for this exposition, flBl ). 

For a generic positive random variable M, let C Ta (X; M) = E[(l + \M)~ Ta ] 
denote its Cauchy-Stieltjes transform of order to. Hereafter we shall refer to this 
as Cauchy transform. Similar to Laplace transforms, C T(T (A;M) uniquely char- 
acterizes the law of M. Let R denote a non-negative random variable with dis- 
tribution function Fr. A random variable M, depending on parameters (t<t, R), 
is said to be a Dirichlet mean of order ra if 



- logC TCT (A; M) = TcrE[log(l + XR)] := Taip R {X) < 00. (2.1) 
Equivalcntly M satisfies the stochastic equation, 

M = f3 TaA M + (1 - j3 Tail )R. 

We denote such a variables as M = M T(J (Fr). Importantly, Cifarelli and Regazz- 
ini |23h (see also apply inversion formula to obtain the distributional for- 

mula for M Ta (Fn). In general these are expressed in terms of Abel- type trans- 
forms. An exception is the case of tct = 1, where the density of Mi(Fr) can be 
expressed as, 

- sin^^e-*^, (2.2) 

7T 

where $r(x) — E[log \x — R\Ifjty x \]. Additionally, 

/•OO 

Wb(A) = to \ (1 - e sA )s- 1 E[e- s /- R ]<is 
Jo 
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is also the Levy exponent of an infinitely divisible random variable with Levy 
density Tcrs _1 E[e _s / fl ]. We say that such a random variable is GGCfYcr, R) and 
may be represented in distribution as a gamma scale mixture 

jraM Ta {F R ) = lT p Ta , T(1 _ a) M Ta {F R ). (2.3) 

Highly relevant to (|2.3[1 . and our exposition, is a result by James (fjjh .that for 
each < a < 1, 

Pr*,T(i-*)M Ta (F R ) = M T (F m J, (2.4) 

where £ CT is a Bernoulli variable with success probability a. Note also that 

A-<t,t(i-ct) = One consequence is that a GGC(tct, i?) variable is also a 

GGC(t, i?^o-) variable. In other words, for a fixed 9 > 0, a GGC(6*, i?) variable is 
a GGC(6>', R^e/9') variable for all 9' > 9. As pointed out in (45), one significant 
point about these multiple representations is that ifO<0 = <7<l, then one 
can set 9' = 1 and use the explicit density formula for Dirichlet means of order 
1, (|2.2p . established by Cifarelli and Regazzini (23) to obtain an explicit repre- 
sentation of the density of such a GGC(er, R) variable. See |45T ) for its precise 
form and further details. 



2. 1 . Identities 

For the case of X ai g, one can show that for 9 > 0, 

C e (X; X afi ) = (1 + A Q )-^ = E[e-™ x <"°], (2.5) 

and for 9 > —a, 

C 1+e (X;X at g) =E[c-^+ lX °-] = (1 + X a )-^ = C e+a (X; X a>e + a ). (2.6) 

We will use (|2.5|) and (|2.6|) to more easily establish the next series of results. 
However, we note that the expressions in (|2.5[) . (|2.6p are not obvious. We will 
provide justification for (|2.5p . when we discuss Linnik variables in the next 
section. Assuming that (12.51) is true, (|2.6[) then follows from an identity due to 
Perman, Pitman and Yor (|75h. 

1 * 0^a t m 

for 9 > —a. (|2.7p is another highly relevant component to our exposition and 
shows that X a , g = l3e +a ,i- a X a j+a- 

Proposition 2.1. The random variables X a ^ are Dirichlet means having the 
following properties: For 9 > 0, 

X a ,8 = Pe,iX a<9 + (1 - /3 8A )X a = Mg(F Xa ), (2.8) 
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and for 9 > ~a and a = (9 + a)/(l + 8), 

X a ,6 = Pe+a,l-aX a .8+ a — P^ e + a l-g ^ce,l+fl = Ml+0 {Fx a ) , (2-9) 

with X a ,g = (3i+e t iX a j + (1 — Pi+e,i)X a ^ a . As special cases of l(2.9\) , 

(U) X a ^i- a = [3l,l-aX a .i = (3^ x- a ^X a ^-q- 

(in) X a = Pa,i-aX ata = /3^™_ Q = M\ (Fx a £ a ) j which yields the identity 

X a = UX a + (1 - C/)X^a 

/or =X a . 

/'too/. follows from Ci +e (A;X Qi e) = C e {\; X a j)Ci{\; X a ). For (JH]), we 

again calculate Ci+g(A; A^g). The first equality is easily checked. For the second 
we use, 

Ci+9 (A; /3^j™ a i-^^Xa^i+e) = Ci+e (A Q , /^e+^i^))- (2-10) 

In order to establish the equivalence to Mi + g(Fx Q £„ ), first note that (|2.8p es- 
tablishes X a , g+a = Mg +a (Fx a )- The result is then concluded by an application 
of (E3D for t = 1 + 9, a = {9 + a) /{I + 9), and R = X a . □ 

The next result establishes results for the larger class of variables defined 
with (12. 3|) and (|2.4|) in mind. Equation (I2.9j) of Proposition 2.1 is an important 
special case. 

Proposition 2.2. Define, for r > and < c < 1, random variables X a a ) = 
f3Ta,T(i-a)X a ,To- Then 

X^aJ = P T a,T(l-a)Xa,Ta = fi^Z T (i- g) ^a,T = M T (Fx a £„ ), 

which leads to the identity 

B 1/a 

Proof. The result is easily checked by following arguments similar to those used 
to establish (|2.9p . Hence we just note that one uses the calculation, C T (A; Xa)-) = 
C T(J (\; X a ^ Ta ), in place of (|2.10p . The equality (|2. 1 1[) follows immediately since 
stable random variables S a are simplifiable, see (|2ll ). □ 

Note that Propositions 2.1 and 2.2 show that 

-logC 9 (A;X a , e ) = -log(l + A a ) = 9E[\og(l + XX a )}. 
a 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. J ames / Lamperti Laws 



S 



2.2. Densities and explicit mixture representations 

We first describe some more pertinent features of X a , see also (|48l ). 

Proposition 2.3. Let X a = S a /S' a , having density Then, 
(i) The cdf of X a can be represented explicitly as 



Fx a (x) = 
(ii) Its inverse is given by 



1 



1 



cot 1 ( cot(7ra) + 



sin(7ra) 



FxM = 



sm(7ra(?/)) 



-i l/a 



sin(7ra(l — y)) 
(Hi) The equations \2.12\l and i2.13\) yield the identity, 



sm(iraF Xa (y)) 



(iv) Additionally, 



cos(naF Xa {y)) 



y a sin(7ra(l - F X(X (y))) 

y a sin(7ra) 
[y 2a + 2y a cos(na) + 1] 1/2 

y a cos(-7ra) + 1 
[y 2a + 2y a cos(?ra) + 1] 1/2 



(2.12) 



(2.13) 



(2.14) 



Proof. This derivation of the cdf is influenced by arguments in Fujita and 
Yor |36l ) where it becomes clear that it is easier to work with the density of 
(X a ) a . Specifically the density of (X a ) a is given by 



sin(7ra) 



1 



7ra y 2 + 2?/cos(7ra) + 1 



for y > 0. 



It it then easy to obtain the form of the cdf of (X a ) a by direct integration. 
Now using the fact that this equates with Fx a (y 1 ^ a ) yields statement [(i)]. 
Statement [(ii)] then follows by using properties of the inverse cotangent. In 
order to establish [(iii)], use (|2.13[) which yields the identity, 



y = F-l(F Xa (y)) 
Hence statement [(ii)] follows. 



sm(ira(F Xa (y))) 



-i l/a 



sin(7ra(l - F Xa (y))) 



(2.15) 



□ 



We now focus on obtaining explicit distributional formulae for the pertinent 
random variables based on their representations as Dirichlet means. In relation 
to this Proposition 2.3 gives precise details on the pertinent cdf F Xa , it then 
remains to obtain a nice expression for the quantity 



§ a (x) :=$ Xa (x)=E[log\x-X a \ 
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for x > 0. The key to calculating $ Q (x) is the fact that we showed that X a 
is a mean functional of the type Mi(F^ a x a ), as described in Proposition 2.1. 
This sets up an equivalence between the form of the density of X a obtained by 
Lamperti (|66T ) and that of Mi(F^ a x a ), obtained from (|2.2|) . Hence we have the 
following calculation. 

Proposition 2.4. For < a < 1, and x > 0, 

$«(x) = -!- log(x 2a + 2x a cos(mr) + 1). 
2a 

Proof. Since X Q = Mi(Fx a £„), it follows by using (|2.2p that the density of X Q 
satishes the equivalence, 

f Xa {x) = -sm(- ! r a [l-F Xa (x)})e- a ^x a - 1 . 

Where on the left hand side we use the expression in (jl.ip . Now applying the 
identity in (|2.14[) shows that, 

fx ( x ) = 1 x^ 1 sin(7ra) c - a$oW _ 

TT [ x 2q + 2a; a cos(Tra) + 1] 1/2 

Solving this expression for $ a (x) concludes the result. □ 
Set 

P a ,r(x a ) = - arctan ( — s * n ( na ) j = ^[j _ p 
a \cos(7ra) + x a / 

and define the function, 

a ( \ xT ~ X sin(p a , T (x")) 

A a>T (x) = :zr . (2-16) 

I" [x 2q + 2x Q cos(a7r) + 1] 2a 

We next obtain density formula for a key class of random variables that 
includes the case of X a , and X a \. 

Proposition 2.5. For < er < 1, and x > 0, the densities of the random 
variables 

with X a = X^l and X a i = X^ \, can be expressed as A a r7 (x), given in \2.16\) . 
Furthermore, 

(i) X ( °l = F x \U a , a ) where U a , a = Fx S ([x£l] a ) has density 



sin(7ra) 



sin(7raw) 



sin(7rcr(l — u)) 
■ ) \, '{,0<u<l. 

Sm(7TQ(l — u) 



(li) If0<a<a, then X^l = [^^-.^X. 
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Proof. The representations of is just a special case of Proposition 2.2. 

The density is calculated based on Proposition 2.4 and the results discussed 
in (IBT ) and (23), as mentioned previously. Statement [(i)] takes advantage of 
the properties of Fx a and is otherwise straightforward to obtain. Statement 
[(h)] is just a manipulation of the beta random variables. □ 

One important aspect of the previous result is that we can use it to obtain 
density/mixture representations for the following Lamperti random variables. 
This is facilitated by the identity (|2.9I ) 

Proposition 2.6. Suppose that < 8 < 1 — a, then a<cr* = 8 + a<l and 
there is the distributional identity, 

y — R y — R v^" 7 *) 

In particular, X a ^ a = 0\,\— a X a ^. 

(i) Hence for < 8 < 1 — a, the density of X a _g, can be written as, 

t Aa '°' iX /Zl du,x>0, (2.17) 
Jo u[l — u) 

where A a ^ a *(x) > 0, is the density of X^' . When 8 = 0, the density is 
A QjCt (a;) equating with hl.l)) . 

(ii) As a special case, when a < 1/2, X a<a = Bi, a X^ 2 *', where X^^ has 
density, 

sin(7ra) 2ax 2a ~ 1 [cos(-7ra) + x a ] 



7T [x 2a + 2x Q cos(aTr) + l] 2 ' 



(2.18) 



Proof. The result follows from Propositions 2.2 and 12.51 by writing 

X a ,9 — @0+a,l-aX a ,6+a = 0l,e0e+a,l—(fi+a)Xa,6+a- 

The simplification in (|2.18p follows from, 

sin(27ra) + 2x a sin(7ra) 



3in(27ra[l - F Xa {x)}) 



1 + 2x a cos(na) + x 2 



□ 



The previous results allow one to obtain simple mixture representations or 
densities for X a fi in the range < 8 < 1 — a, and 8 = 1. The fact that we 
obtain such results for a continuous range of 8 is significant, as shown in the 
next result. 

Proposition 2.7. Set 8 = Ylj=i ®j where 8j > 0. Furthermore, let (Di, . . . , D^) 
denote a Dirichlet random vector having density proportional to Yii=i x i* ■ That 
is each Di = Then, 



X a ,6 — DjX ct: g j 



3=1 
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where X a ,g j are mutually independent and independent of (-Di, . . . , Dk). When 
9j are chosen such that < 6j < 1 — a, each X aj g j has an explicit density fx a a 
described in \2.11^ . When 9 = k, one can use 9j = 1. 

Proof. Since we have shown that X a ^ = Mg(Fx a ), this result follows directly 
as a special case of Hjort and Ongaro( (41 ), Proposition 9). □ 



2.3. Positive Linnik variables 

For 9 > 0, 

X a ,e=lTs a (2.19) 

denotes the class of generalized Linnik variables as considered in {7(1 EH; 
281 : 68t 43). The results in the previous section depend on the validity of the 



transforms in (|2.5p and (|2.6p . It is evident, and known, that the expressions is 
the Laplace transform of Xa,e for 9 > 0. Hence (|2.5I) is verified if one shows that 

that Xa,e = leX a fi, for 9 > 0. It is already known, using a result of Devroye(27) 
combined with (|2.7[) . that 

Xa,a = l\ S a = "flX a = J a X a , a . (2.20) 

Furthermore, from Bondesson ((|16l). p. 38), it follows that Xa,e are GGC(#, Fx a ). 
In the next result we will verify the usage of (|2.5p . and use the X^\ to obtain 
explicit density representations. In this regard, it is is important to note that 
we do not need explicit results for X a ,g to get corresponding results for Xafi- I n 
addition we obtain some interesting identities. 

Proposition 2.8. For all 9 > 0, \a,e is a GGC(9, Fx a ) variable that satisfies, 

Xa.fi = leX a fi. 

For < 9 = a < 1, Xa .M = 1\X£{. and hence has the density, 

poo 

fx.A*)= / ^ ,y y- l ^ a Av)dy- (2.21) 

Jo 

See also \2.2J$ for 9 > 0. Additionally, 
(i) For 9 > -a, 

Xa,9+a = Je+aXa.g+a = ~fi + gX a ,g. (2.22) 

(ii) Hence, for 9 > —a, 

1/q _ 16+a 71+9 / O o\ 

(Hi) For -a < 9 < k, k = 0,1,2,..., 



Xafi+a — lk+l^a,kP^e+a k- 
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(iv) For 9 = Yli=i ®i > Oj Xq,0 = SiLiXcMi? where Xa.Si are independent. 
Proof. From ()2 . 1 9[) and using the identity 

e~£ =E[c-? s °], 
it is easy to see that the density can be expressed as 

/*°° a 







oc 



oc x 9 - 1 / / e-^{v/s) e v- e f a (v)f a (s)dvds 
Jo Jo 

yielding the equivalence with "fgX a ,g. The expression in (I2.2ip is due to Proposi- 
tion 2.5. Statement [(i)] follows from (|2.7p . Statement [(ii)] follows by removing 
S a , which is justified since it is a simplifiable variable. For (iii), apply Proposi- 
tion 2.2. [(iv)] follows from infinite divisibility. □ 

Remark 2.1. It is not difficult to show that a general expression for the density 
of Xa.9, f or oil 9 > 0, is obtained by replacing A Q O - by A Q .# as follows, 

/x« 9 (*) = -f ^^ 0F ^ d \, (2.24) 
7T J [ x 2a + 2x a cos(air) + 1] 2 ° 

However, A aj g can take negative values when 9 > 1, so this does not in general 
yield a mixture representation for Xa.e- Nonetheless, it may not be difficult to 
evaluate numerically, which is relevant for the Mittag-Leffier functions discussed 
in section 3. 

Remark 2.2. The gamma identity in statement i2.23\) of the previous proposi- 
tion is quite remarkable, and as we shall see below has some interesting implica- 
tions. We note that although not obvious, our result coincides with a variation of 
Bertoin and Yor (XlA). Lemma 6). Checking moments one can see that, in their 
notation, J SyS/a = S~% for s > 0. and for 9 > -a, Ji+g^ g+ / a y a = SlS. Our 
work provides some additional interpretation of their variables, see also See 
Kotlarski Xdh ) for a general characterization of cases where products of variables 
result in gamma variables. 



2-4- Generalized Pareto Laws 

Influenced in part by (|2.23p , we next look at relationships between the Lamperti 
laws and a class of generalized Pareto distributions. We note that the next 
result also plays an important role in section 4, when discussing continuous 
state branching processes. Define random variables, 

wV<* ( Uf \ " 
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These represent a sub-class of generalized Pareto distributions with cdf and 
density given as 

y e Oy 8 - 1 



w (i + 2/ «)«/«' J <: e vy; (i + y «) (9+Q)/Q ' 

Proposition 2.9. LetU denote a Uniform[0, 1] random variable, then for 9 > 0, 
(i) there is the identity, 



l-Uv J V 7i / 7i 7i 

(ii) For < er < 1, i/ie random variable S QjCT = 71 /I^'j, /ias Laplace trans- 
form, 

E[ e - AS "-] = 1 - A ff (l + A Q r CTA \ (2.25) 
Proof. Statement [(i)] is an application of (|2.23[) For statement [(ii)] notice that 

P(-^_ > A) = E[e- AS °-"], 
but this is the survival function of the random variable 

7l y (a) £ la S a £ i /a 
71 b a . a 71 

□ 



Remark 2.3. j4s we shall discuss in section 4, Statement [(H)], 12.251) serves 
to identify explicitly the (unknown) limiting distribution obtained by (92) and 
ftv\ ) corresponding to a = 1. It is relevant also to note that the density of 
W^/i is the only case that corresponds to a Laplace transform. So here we see 
a distinguishing feature of X a \. 

2.5. z -variables and hyperbolic laws 

Proposition 2.9, along with the works of |l4; 84; 95; 4^), motivates us to con- 
sider several questions related to z-distributions, which are distributed as the 
logarithm of independent gamma variables. We also believe that some of the 
variables we discuss will be of interest in terms of applications along the lines 
discussed in Q) and Huj42l). In fact, (|42T ) suggests the use of a class of variables 
that turn out to be equivalent in distribution to log(Jf a ). Naturally we do this 
in the spirit of highlighting what one can do with Lamperti laws. We also obtain 
additional information about these variables. 

In particular, for illustration, we consider the following generic type of prob- 
lem. Suppose for generic variables X, Y, Z with Z and Y independent there is 
the following relation 

X = Y + Z. 
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One natural question is to ask given explicit information about X and Y , what 
Z satisfies the above equation? In addition does Z have an explicit density or 
mixture representation? Notice also that if the characteristic function of Z is 
not known then we can use X and Y to obtain this. We will consider Z that 
are variants of Lamperti laws. 

We first give a brief discussion on z-distributions and cite some relations 
to log(Xa) that exist in the literature but are probably not well known. Fol- 
lowing (|95l ). the class of z-distributions are defined as 7r _1 log(7e 1 /7e 2 ), having 
characteristic function 



r(0i)r(0 2 ) 

As special cases, the variables, for < a < 1, M a = tt^ 1 log(7 CT /7i_ CT ), have a 
Meixner distributions with characteristic function 

E[e"^] = fff'H (2.26) 
cosh(A — is a ) 

where e a — ir (a — 1/2). Note that a Meixner distributed is usually defined as 

(1/2)M CT . §i = 7T _1 log(t//(l — U) has a logistic distribution with characteristic 
function, 

E [ e iASil A 



sinh(A) ' 



and Ci = 7r 1 log(7i/2/7i/2) nas the hyperbolic distribution with characteristic 
function. 

1 



r e ttCn _ 



cosh(A) 

It is known that the variables Si and Ci satisfy 

Ci=Si+Ti, (2.27) 
where Ti is an independent variable having characteristic function 



E[< 



iXT _ tanh(A) 



' A 

Biane and Yor flil ) showed that the density of Ti is 

1 4 
/ti(^) = — log(coth(— \x )), — oo < x < oo. 

7T 7T 

In regards to this and other considerations note that the characteristic func- 
tion of cot" 1 log(X a ) is equivalent to 

E [ e ^iog(x Q)] = sinh(aA) _ 
a sinh(A) 
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This expression can be found in Chaumont and Yor ((|2lh. p. 147). In addition we 
see that this characteristic function agrees with the class of generalized secant 
hyperbolic distributions discussed for instance in (I42|). It follows from the first 
Laplace transform given in Biane, Pitman and Yor f(|14T). p. 459, eq. (6.1)), that 

wn^ 1 \og(X a ) = B(T"(R 3 )), where B{-) is a Brownian motion and, independent 
of £>(•), T 1 Q (i?3) is the first hitting time of 1 by a Bessel process i?3 of dimension 

3 starting at a. From Chaumont and Yor (j2ll . p. 155 and p. 169), lim Q ^o = 
U/(l-U). Hence, 

lim an- 1 log{X a ) = lim B{T?{R 3 )) = B(T?(R 3 )) = S x . 

Note however that tt" 1 log{X 1/2 ) = d = B(T?(Ri)), for i?i a Bessel process of 
dimension 1. 

Now using (|2.7p and (|2.1ip , one can write 

log(X a ) = log(^)+log(§^) (2.28) 



= M^) + -log(^) 



D a,l 



l0g(Y Q , 1 ) + -l0g(/3 (1 i^a). 



Notice that when a = 1/2, §i = 7r 1 log(S' 1 /2,i/2/*S'i/2 1/2)' hence the first 
equality in (|2.28|) coincides with (I2.2T[) suggesting that this line is the natural 
generalization. Furthermore, it follows that 

Tl = 7T l g(- ). 



/2.1/2 



This illustration can obviously be applied to other variables using (12 . T[) and 
(|2.11[) . However we will not attempt to develop a thorough analysis here. Rather 
we will present two results that give a sense of what else can be done. In the first 
result below we are able to obtain a description of the characteristic function of 
X afi andX<g. 

Proposition 2.10. From Provosition \2.9l it follows that, 
(i) for 9>0, 

— log f — ^ = log f — ^ +log(*„, e ). 



(ii) For 9 > —a 



r(^+*^)r(l~^)r(l + fl) 

r(i + + iA)r(i-£)r(*£s) 
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- log f^u log r^u ^(41) 



(Hi) For < a < 1, 



a \7i/ \7i 
where log(X^ CT ]), has density 

1 sin{p a , a (e za )) 



7T [ e -2 2Q + 2e~ za cos(a7r) + 1] 2 ° 
and characteristic function 



-co < z < oo. 



sinh(A^) r(f + g)r(i - f) 
a- r(*) 



Proof. This follows as a simple consequence of our previous results and the 
characteristic functions of z-distributions. □ 

The next result illustrates some cases where the characteristic functions are 
based on hyperbolic functions and all random variables have explicit densities. 
Define, 

tt A a-i A ^g-.i-g- v 

° - -JJ^- - a 

P\-a,a Pl-a,(T 

where the equality follows from (|2.11ll . In addition for aS < 6 < a(l — S), for 
5 < 1/2, define 

1_ 

(6) A ( \ " 5q,i-e 



where one can easily check that the density of S a ,i-e/S a! e, denoted as fi-e.e, 
satisfies 

fi-e,e( x ) = c a,ex~^~ 6 ' fxc^ix) = c a fix~ {1 ~ 8) A a ,i(a;), 

for 

_ r(i/a + i)r(g + i)r(2-g) 
Ca - e ~ r(0/a + i)r((i - 0)/a + 1) ' 

Proposition 2.11. For < a < 1, 5 < 1/2 and a<5 < 9 < a(l — 5), f/iere are 
t/ie following relationships: 



(i) ilog(^) =log(^) +log(ff 0)<7 ). ffence, 



gr j«> losCiJc »)i cos(e CT ) sinh(Aa) sinh(aA) cos(eo-) sinh(A) 
Aa cosh(A — ie a ) asinh(A) Acosh(A — ie a ) 

where e„ = Tr(a — 1/2). 
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(ii) Ilog(^) ilog(^) +log(4 5) e ). With, 

jr^iogii™,), = cos(ea) cosh(Aa - zee) 
cos(ee) cosh(A — ies) 



(Hi) When S = 1/2, then 6 = a/2, and L^£j 2 = S al _ a / 2 / 'S aa / 2 . wh 
lo &( L a,a/2)> haS densit V 

c a ,a/2 e z % sin{p aS (e za )) 
■ ■ j-, — oo < z < oo, 

TT [ e 2za + 2 e za cos ( a7r ) + 1] 5S 

and characteristic function 

Efe ^ togfi^), = cosh(A - ie a/2 ) 



ere 



cos(e a / 2 ) cosh(A/a) 

Proof. All the characteristic functions follow from that of Si and (|2.26p . In order 
to establish [(i)] we use the identity ji/S a — ^ a and apply this to the second 
definition of H a cr . Statement [(ii)] follows from a manipulation of (|2 . 23[) to force 
the form of the first two variables. The choice of 1 — 9 and 9 in the definition 
of was deliberately made so that we could get an explicit expression of the 
density of the relevant ratios of stable variables. □ 



3. Mittag-Leffler functions 

In this section we obtain integral representations, and other identities for a 
generalization of the Mittag-Leffler function given by 

E^(-A) = Y ( ~y* + for 9 > -a (3.1) 

Q < 1+e ^ Q k\ T(ak + 9 + l) 1 ; 

where 

T(9/a + l + k) 
(6/a + l) k = v{e/a + 1) • 

So when 9 = 0, one recovers the Mittag-Leffler function as, 

E Q , 1 (-A) = Ei, 1 (-A) = Ei 0) (-A). 

The function (|3.ip is a special case of the function introduced by Prabhakar jHI), 



w-t^rgh (3 - 2) 

where (p, ^,7 £ C, Re( i o) > 0). That is the case where 7 = (9 + a) /a and 
fi = 9 + 1. The quantity (|3.ip represents a special sub-class of yet more general 
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Mittag-Leffler type functions which are discussed, for, instance in Kilbas, Saigo 
and Megumi See also (@; 0; H; 0; S H H [H H2) . 

Note that using simple cancelations involving gamma functions it is easy to 
show that for 9 > 0, 



? {e/a+i) ( _^ _ r(6>) (a/a), 
r(l + 6>) 



ELT + y } (--)-Tv^Cr(--)- (3-3) 



Recall that the Mittag-Leffler function can be expressed as the Laplace trans- 
form of S~ a . One can show that by taking a Taylor expansion and calculating 
moments that 

E[ e ->/sS, 9] = T{0 + l)E^t 1} (-^ (3.4) 

One consequence of this observation is that one can use a Monte-Carlo method 
based on S a ^g to evaluate this quantity. The next result develops more connec- 
tions with X a fi and Xa,e- 

Proposition 3.1. For 9 > —a, 

(i) fl/X ae (x) = X- fx a , e {x) 

fti) fy 1 /x a Jx)=T(9 + l)x- f Xae+a (x) 
(Hi) For 9 > 0, 

„\ jrp\ —X.X a qI T?{®f a 



>x)= E[e-* x «,°] = K' a \-x a ) = mx^f^ 



Hence if 9 = a then these expressions are explicitly determined by \2.21\) . 
otherwise one might use {2.2$ - 
(v) Applying Proposition 2.7, it follows that for X)j=i @i = ®i where 9i > 0, 

where S k = {(x 1: . . . ,x k ) ■ < x i < !}• 

The result is fairly straightforward using the basic definitions as ratios of 
stable variables and the identities in Proposition 2.5. We omit the details. 

Remark 3.1. Recall that for a = 1/2, Xi/2,0 = 76+1/2/71/2- Using the notation 
in Pitman ((7§), eq. (88), eq(98) and Lemma 15), and noting ff79l), eq. (104)), 
the conditional moments of the meander length (1 — G\), of Brownian motion on 
[0, 1] conditioned on its local time L\ is related to the generalized Mittag-Leffler 
function as follows, 



5i| e+1/2 ]/M20+i)(A) 



E 1 / 2j o[(l-Gi) 0+1/2 |ii = >/2A] = E[e~* 

9+1/2] , 

-(20+1 

= n(0+l/2\\\), 
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which corresponds to the moments of the structural distribution of the Brownian 
excursion partition. h-2q(X) is a Hermite function, here q = 9 + 1/2. 



4. The explicit Levy density of Stable CSBPs and the 
Zolotarev-Slack distribution 

We now show how our results lead to am explicit identification of the Levy 
density of the semigroup of stable continuous state branching processes of index 
1 < S < 2, that is 8 = 1 + a and the limiting distributions first obtained by 



Zolotarev(|97T) and Slack (|87l ) . We also mention briefly its connection to the work 
of H; 0) on beta coalescents. From Lamperti, (|o5l ) continuous state branching 
processes (CSBP) are Markov processes that can be characterized as limits of 
Galton- Watson branching processes when the population size grows to infinity. 
A stable (CSBP) process (Y t ,t > 0) is a Markov process whose semigroup 

is specified by, 

E a [e" AYt ] = E[c- XYt \Y =a] = e~ av *W (4.1) 

where 

/>OG 

vt= / (1 - c- sX )v t (ds) = X(at + X a y 1/a . 



Furthermore, see for instance (0;@), there exists a process (Y(t, a); t > 0, a > 0) 
such that for each t, Y{t, •) is a compound Poisson process with intensity v t . 
Related to this result, Kawazu and Watanabe (|53l ) show that all continuous 
state branching processes with immigration arise as limits of Galton- Watson 
processes with immigration. Analogous to (|4.1|) . Theorem 2.3 of their work yields 
the limiting (1 + a) stable (CSBP) with immigration (Y(t), t > 0) satisfying, 

E a [e-^} = (1 + a\ a t)-^E a [e- Yt }. (4.2) 



It is evident from (|4.2p , that the entrance laws of Y are positive Linnik distri- 
butions. However, the intensity v t , which plays a fundamental role in (0; 0) is 
only known up to its Laplace transform. It is known that this Laplace transform 
coincides, up to some scaling factors, with the Laplace transforms of the limit- 



ing distributions obtained by Zolotarev (|97l ) and Slack (|87l ). Before identifying 



these expressions we recount the limiting distributions obtained by Slack 
and Berestycki, Berestycki and Schweinsberg {!). As described, for instance, 
in (Q), Slack's result describes the limiting distribution, say fj, a , of the number 
of offspring in generation n of a critical Galton Watson process, rescaled to have 
mean 1 and conditioned to be positive, when the offspring distribution is in the 
domain of attraction of a stable law of index 1 < S < 2. This result complements 



Yaglom's (|92l ) well known result for the case where the offspring distribution has 



finite variance. In that case the limitingdistribution is exponential with mean 1 



Precisely, following the exposition in (|72[). we state a variation of Slack's result 



Proposition 4.1. (Slack(1968) (87)). Let Z = (Z n ,n > 0) denote a super- 
critical Galton Watson process initiated by a single process. Furthermore, sup- 
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pose the non- extinction probability Q n = P(Z n > 0), satisfies 

where L(x) is a slowly varying function. Then, 

lim P(Q n Z n < x\Z n > 0) = (i a ({0, x}) (4.3) 



n — >oc 



where for each < a < 1, fJ, a is the distribution of a random variable S Q 
satisfying 

e- Xw n a {dw) = E[e- AS «] = 1 - A(l + \ a y 1/a . (4.4) 

Zolotarev f (|97l ). Theorem 7) also obtained this limit in the case of a class of 
continuous parameter regular branching processes. However, prior to our work, 
an explicit description of its density or corresponding random variable was not 
known. It is now evident from (|4.4|) that S Q = in Proposition 2.9. as we 



mentioned previously, These limits and or discussions related to (|4.1ll . (|4.2| 
appear more recently in for instance, (6; 9; 33; 63; 64; 74). Before we summarize 
our results we shall say a bit more about the context of (@). Random variables 
with law _u a arise in the work of Berestycki, Berestycki and Schweinsberg (@), 
see also(|9|), in connection with Beta (2 — 5,5) coalescents for 1 < 5 < 2. Sec 
in particular ((0), Theorem 1.2.) Equivalently these are Beta (1 — a, 1 + a) 
coalescents. 

In addition, there is a related result of Berestycki, Berestycki and Schweins- 
berg (0) that, with some work, would have otherwise allowed us describe the 
law \i a . We quote their result below, 

Proposition 4.2 (Berestycki, Berestycki, and Schweinsberg ((0), Proposition 
1.5)). Let (n(i), t > 0) denote a Beta (1 — a, 1 + a) coalescent where < a < 1, 
and let K(t) denote the asymptotic frequency of the block o/IT(t) containing 1. 
Then 

{T(a + 2)t- 1 )"K{t) -> Ca as 1 10, (4.5) 

d 

where £ a is a random variable satisfying, 

E[e- A H = (l + \ a )- (1+a}/a . (4.6) 
Furthermore, as noted in (0), £ a has the size biased distribution 

F(C a edx)=x(jL a (dx). (4.7) 



We now summarize our result which again demonstrates the relevance of X a ^\ 

Proposition 4.3. 1 

density A a i. Then, 



Proposition 4.3. For < a < 1, recall that X a i = S a /S a ,i, has explicit 
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(i) £ Q , the random variable described in |^.5[ j and \4-@i) , satisfies 

(a = ll^-a.l — Xa,l+a- 

(ii) Let Yi a and /i a be as in Ij4.3\ l and fi4-4\ )> then 



„ ,4 71 d_ v 

fmj Furthermore, for each x > 0, 

P(S Q > z) = A* a ([s,oo)) = E£/ a) (-z Q ) = / 7l x a>1 (z). 

(M/J Sei # t = , t/ien £/ie Lew?/ density v t corresponding to the (1 + a) 

(CSBP) specified by 0771) is, 



(at) ^"V- 1 ( i±g) x", 
" t(a;) = r(l + a) E "- f +" ( '^' 



= (ai) - / e C*^yA a ,i(t/)dy. 
Jo 

Proof. Statements [(i)] and [(ii)] are now quite obvious from Propositions 2.8 
and 2.9. Statements [(iii)] and [(iv)] are deduced from Propositions 2.5 and 3.1, 
along with the specifications for v t , in terms of 9 t described in ((@), Lemma 
2.2). □ 

5. Occupation times of generalized Bessel bridges 

We now show how our results for X a j can be used to obtain new results related 
to At~ , which is equivalent in distribution to P a .e(p) : under P^ P g. This can be 
seen as a continuation of a subset of the work of (|47T ). who looked at more 
general PD(a,6 l ) mean functionals , where, with the exception of a — 1/2, the 
best results for describing the density of P a fi{p) were obtained for 9 = 1, and 
9 = 1 — a. The results for a = 1/2 are classic. For a — 1/2, and p = 1/2, 
Levy ([63) showed that Af under (1/2, 0) and (1/2, 1/2), follow the Arcsine and 
Uniform[0, 1] distributions respectively. A general formula for (1/2,0), for all 
9 > —1/2 can be found in Carlton (p~9h . see also |H3). We also obtain results for 
time spent positive on certain random subsets of [0, 1], and also develop some 
interesting stochastic equations. As a highlight, we obtain explicit results for the 
case of 6 = a, corresponding to the the time spent positive of a Bessel bridge 
, 1 1 . In this case the best previous expressions were obtained independently 



on 



in (|47|;|93D 



Consider now the following stochastic equations and Cauchy transforms that 
can be found in |47h with further references, for 9 > 

P a ,e(p) = 0e,iP a ,e(p) + 0-- Pe,i)P a fl{p) (5.1) 
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and for 9 > —a, 

Potfiip) = Pe+a,l- a P a ,e+a(p) + (1 - 00+at,l-a)tip- (5-2) 

Additionally there are the Cauchy transforms for 9 > 0, 

Ce(X; P a ,e{p)) = (?+(! + A)"?) - * - e - fl *&W (5.3) 

where 

^(A) = E[log(l + \P a , (p))] = E< rt [log(l + XAf)}, 
and for 9 > —a, 

Ci+e(A;P Q , e (p)) = 5+^— (5.4) 

(«+(l + A) a p) — 

= C6»+a(A; P a ,8+a{p))Cl-a(X; £p). 

The first equation fl57Q) shows that P Q ,e(p) = M e iF P n(rA ) for > 0. The 
second equation (|5.2p , see for instancelfioh. |29l). (44), (|80h. for some other inter- 
pretations and applications, can be traced to Pitman and Yor ((|8l|), Theorem 
1.3.1), and Perman, Pitman and Yor ((f75h . Theorem 3.8, Lemma 3.11) as fol- 
lows; Let Aq i = J 1 1(s s >o)^ s denote the time spent positive of B up till time 

Cm, which is the time of the last zero of B before time 1. Then under P^ p l there 
is the equivalence 

Pa,e+a(p), Pe+a,l-a) 

This shows that (|5.2p can be rewritten in terms of the following decomposition, 
A+ £ Aq i + (1 - Gi% = GtA^ + (1 - Gi)f p . 



See for example Enriquez, Lucas, and Simenhaus |31j) for an interesting recent 
application of this expression. 

We now show that the density of P a ,e{p) can be expressed in terms of the 
density of X a ^. Hereafter, define r p (y) = y/(c(l — y)), for c a = p/(l — p). 

Proposition 5.1. For 9 > —a, let R a ^g — cX aj g/(cX a ^ + 1). Then 

^](A+ G dy) = {1 ~ y)9 P(R a , g G dy) 
(l-p)~ 

where A^ = P a ^(p). Hence as special cases, using Propositions 2.5 and 2.6, 

(1) Pjl(A+ G dy)/dy = (1 - y)- x p- l /<*A aA (r p (y)) = fl^y). 
(ii) For < 9 < 1 - a, and a* = 9 + a, 



% 6 {At G dy)/dy = ^ x / 

where A Qi(T »(x) > is i/ie density of X^ . 
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Proof. From (|1.3p it follows that for measureable functions g, 
E[g(cX a , e )} = ^^ E%[g(^){A^)- s 



1 T i 



The result is concluded by showing that 

is equal to (1 - P y 0/a E[g(R a j)(l + X afi )~ 6 ]. But this follows from n = (A+/Ai~ ) + 



Pitman and Yor ((|83h. Proposition 15) establishes an interesting relationship 
between the densities of Af and Aq under the law Pffp. Making no changes 
to the essence of their clever argument, one can easily extend this result to all 
(a,ff). Combining this with Proposition 5.1 yields the relationships, 

^ e {A + Gi edy) = ±-?Lp^ g (Aiedy) (5.5) 

(1 — i ) 1+e 
= ~ ^r^F(i? Q . e e dy). 

(1-p) — 

Recall that under Pjfl, Aq = f3e+a,i- a Pa,8+a{p)- The next result describes 
interesting properties of generalizations of this variable. 

Proposition 5.2. For r > and < a < 1, let V ^ ft,^ r{i- a ) < and hence is 
a Dirichlet mean satisfying the stochastic equation, 

V = Pr/ a ,lV + (1 ~ r/a ,l)&, 



d 



for V = V. Then for < p < I, 

Pt„,t{i-*)P<x,t«(p) = P a ApV) = M T (F Pai0(pM ). (5.6) 
This leads to the stochastic equations, for < p < 1 

P a ,ApV) = &,iP' a ,T(P v ') + i 1 ~ Pt,i)P*,o(p&) (5-7) 

= $r,lK >T (pV) + (1 - 0r,l)Pa,o(j>V). 

In the first expression P' a T (pV) denotes a random variable equivalent only in 
distribution to P QjT (pV^). However in the second equation V is the same variable. 

Proof. First note that fi TrJ ,T(i-a)Pa,Ta{p) = M T (F Pa ( p ^)), follows from (|2.4p . 
Note also by the definition of P a ,o(p) it is easy to see that P a ,o(p)£,<r = Pa,o(p£a)- 
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In order to establish the rest of (|5.6|) we can check Cauchy transforms of order 
t using (|5.3[) . For the variable appearing on the left of (|5.6[) this is easy to 
calculate. Applying this to P a ^ T (pV) conditioned on V, its final evaluation rests 
on the simple equality 

(1 -pV + (1 + X) a pV) = 1 + [(1 + X) a - 1} P V. 

Taking the Cauchy transform of order t / a for V yields the result. The second 
equality in (|5.8|) is then due to (|5.1[) . □ 

Next, is one of our main distributional results, which is an analogue of Propo- 
sition 2.6, but also highlights the role of various randomly skewed processes. 

Proposition 5.3. For < a < 1, set R { °\ = cX£l/(cX£l + 1), then the 
density of 

Pa,{l-a)Pa,a{p) = Pg.l (pP { ^ (1-*) ) ) 

is for < y < 1, equivalent to (1 — p) a (1 — y)P(R a [ £ dy)/dy, and is given 
explicitly as, 

^ [y 2 V + 2gpy«(l - y) Q cos(aTr) + (1 - y) 2 V] 2 ° 

In particular f2 Q ,i(y) is £/ie density of P a ,iip) an d ^aAv) * s ^ e density of 
0a,l-aPa,a(p)- In addition, 

(i) if Q < 9 < 1 — a, then for a* = 9 + a, 

0e+a,l-aPa,6+a(p) = Pn^+f)(v0{o* i-<*}) = 01.f)Pni.l(p0^* ^ i-<r* ))i 

and, using II5.6\) , there is the explicit formula determining the densities of 
Af and , 

P i a P l(A+ i edy)/dy = l^F^(A+edy)/dy (5.9) 

1 On a , a .(y/u), 
du. 



o u(l — u) 1 

When 9 = 0, 15. 9j) is Cl a , a {y) w hich agrees with Pitman and Yor((8S), 
Proposition 15). 

(ii) As in Proposition \2. 7[ for any 9 > 0, set 9 — 5Z,- =1 0j /or some integer k 
and 9j > 0. TTiis /eads to the representation, 

k 

p a Ap) = Y, D j p *fM 

3 = 1 

for independent variables P a ,8j(p) and (D\, . . . ,Dk) a Dirichlet vector as 
in Proposition \2. 7| When 9j are chosen such that < 9j < 1 — a each 
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variable has ¥^ g (A^ G dy) given by \5. 9\) with cr* = 6j +a. It suffices to 
choose 6i = 9/k, for < 9 < k(l — a). If 9 — k, one may set 9j = 1 and 
use 

Proof. The various representations of the random variables are due to Proposi- 
tion 5.2 and otherwise an application of the beta/gamma calculus. The density 
Oq, )(t is obtained from A QjCr , which is justified by the exponential tilting rela- 



tionships discussed in James f (|45l ). section 3), see also (|46f ). □ 



Remark 5.1. The random variable P a , T (pV) described in Provosition 1 5. 2\ has 
law 



P(P a ApV) G dx) = P%p(A+ G dx) := / P^(A+ G dx)h 

Jo 



(u)du. 



That is, it may be read as the time spent positive up till one of a process B 
whose excursion lengths, conditional on V, follow a PD{a,r) distribution and is 
otherwise randomly skewed by pV. See also Aldous and Pitman (IA), section 5.1) 
for connections with T-partitions. This is made clear, as follows; For (L t ;0 < 
t < 1) governed by PD{a,9), and letting L t — L t /Li, there is the equivalence, 



P a ,e(u) = inf{i : L t > u}, < u 



< 1. 



In other words, letting P^ e (•) denote the random quantile function of P a ,e, it 

follows that Lt = P^ ^ (t),0 < t < 1. See the next section, section 6, for more 
general V. 

Remark 5.2. We note that in reference to Propositions \5.2\ and \5.3[ set- 
ting Q atT (a,p) — fJ T<T ,T(i-a)Pa,ra(p) leads to a well defined bivariate process 
{Qu,t{o~,p) ■ < p < 1, < a < 1), that has some natural connections to the co- 
agulation operations discussed in Pitman f7c\ ). This observation may be deduced 
from the subordinator representation given in Pitman and Yor ((82). Proposi- 
tion 21). When p = 1, Q a . T (o~, 1) is a Dirichlet process, which corresponds to 
the operation of coagulating PD(a,r) by PD(0, r/a). In general one may write, 

Q a , T {cr,p) = Po, T (cr)P a ,Tcr(p) = P a , T (pP 0l z.(o-)). 

We will not elaborate on this here except to note that in connection with results 
for the standard U-coalescent, setting t = a, p = 1 one recovers (%7&), Corol- 
lary 33, and Proposition 32). In a distributional sense, that is without the nice 
interpretation, one also recovers (f7c\). Corollary 16) by setting a — a = e _t 
and t = l,p = 1. When p ^ I, we expect that one can obtain new, but related, 
interpretations of Q atT (a,p). 
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5.1. Some special cases 

Note that, as in |47l ) (combined with Proposition 5.2), one can rewrite (|5.2[) as, 

a ,1 —a 

= f > a,l+e(p/3(9±<S ) i^))(l - £p) + (1 - fa,l+e(9/3(«±2,i=H)))ff>- 



Besides giving an alternate mixture representation in terms of easily interpreted 
random variables, this also suggests that one can obtain the density of P a j(p) 



if one knows the density of P a ^e+ a (p). In (|47l). it was noted that this could be 
applied for 6 + a = 1, which yields an expression for the density of P a .\- a (p). 
In view of Proposition 15.31 we see that such a density representation can be 
extended to any < 9 < 1 — a. Of course in terms of a density representation 
this is not as good as the expression one can obtain from (|5.9p . since it would 
have to be used twice. In this section we look at some specific cases of random 
variables that have cither appeared in the literature or we anticipate might be 
of some interest. 

Example 5.1.1 (fa, 1 — a) A distribution relevant to phylogenetic models). 



As noted in (|47l ) the case of P Qi i_ Q (p) equates in distribution to the limit 



of a phylogenetic tree model appearing in ( (|38l ). Proposition 20). Here using 



Proposition 15. 31 we obtain a slight improvement over the density given in [ (|47l ). 
Corollary (6.1)]. Since under ¥^\_ a , A^, = (3i^ a P a _i(p), we have, 



6 dy)/dy = ^—^P^ a (Atedy)/dy (5.10) 
1 (l-a)n a ,(|//«) 



1 



o 



it(l — u) a 



-du. 



Example 5.1.2 (The case of (3i +a ^- a P a ^ +a (pj). 
Under P&J. 
is given by 



Under Pj^, = l3i +a ,i- a P a ,i+a{p) = P n ^(p0^+ a ^-"))- Hence its density 



^\{A+ i edy)/dy = ^-^ ) £l aA {y). 

In view of the literature related to section 4 we believe this variable will be of 
interest. 

We now address some harder cases. 

Example 5.1.3 ((a, a) Occupation time of a Bessel Bridge). 
Obtaining density expressions for the general case of A^ when B is a Bessel 
bridge has been difficult, except for the case of a — 1/2. Due to the importance 
of the PD(a, a) family this quantity arises in many contexts, see for instance Al- 
dous and Pitman |l|). The best results were obtained independently by Yano(93) 



and (|47l). who give expressions in terms of Abel type transforms. That is to say 
integrals of possibly non- negative functions. Hence this does not yield a mixture 
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representation for P a>a (p) = A*. Here we show how our results in the previous 
section can be used to achieve this. Under P$ a , 

j4(3i = 02a, 1- a Pa, 2a (p) = Pg 1 +a (pP(2, ) )' 

Hence for a < 1/2, we can apply statement [(i)] of Proposition 5.3 writing, 

A% 1 = Pl, a P a ,l(pP(2,l=£i)) 



to get 



Where 

^a,2a(j/) 



»• dy)/dy = l—Lv^(A+ e dy)/dy (5.11) 

1 aVL a ^ a {y/u) 
o «(1 — u) 



2sin(7ra) py 2a 1 {l — y) a [qy a + cos(7ra)p(l — y) a ] 
7T [ 2/2Qg2 + 2 gro «(l - y) Q cos(Tra) + (1 - yfVf 



is the density of P a ^(pf3^ 2 i-^ O- 

When a > 1/2, we, at present, need to resort to statement [(h)] of Proposition 
5.3. So, for instance, for a < 3/2, it follows that 

Pa,a(p) = f3 a /2,a/2P a ,a/2(p) + (1 - P a /2, a /2)P' a , a /2^P) 

where P a , a /2(p) an( l P' a a /2^P) are ^ variables having distribution JP^m £ 
dt/) obtainable from (|5.9p . 



Example 5.1.4 ((a, a — 1), and fragmentation equations.) 

Suppose that we are interested in the case of 9 = a — 1, that under P^; 

Pa,a-i(p) = -^i*- Of course this only makes sense for a > 1/2. Notice that 



Aq 1 — 02a-l,l-aPa,2a-l(p) 

so we can apply (|5.9p directly if a < 2/3. It is then interesting to note what 
other quantities we can obtain. We can use another stochastic equation that 
takes the form, 

Pafiip) = Pe+aS,l-asPa,9+a8(p) + (1 — Po 

for 9 > —ad, < 6 < 1. Note that this equation is not well known but it is 
simple to check. Furthermore, a close inspection shows that it is a nice way to 
code Pitman's f78|) fragmentation. As as special case, set 6 — (1 — a) /a and 
9 = a, so we obtain 

P a ,a(p) = Pl,aPa,l(p) + (1 - Pl,a)Pa,a-l(p)- 
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Another consideration is the choice of 9 = 1 — a, that might be relevant 
to dH), yields 

P a ,l-a(p) = fh(l-a),aPa,2(l-a)(p) + 0- ~ P2(l- a ),a)Pa,a-l{p)- 

6. Power scaling property and randomly skewed processes 

We saw that in the previous section the random processes P a . T (pV), where V is 
a beta variable occurs naturally and plays an interesting role. An interpretation 
in terms of occupation times of randomly skewed processes is mentioned in 
Remark 5.1., and an interpretation via coagulation processes is hinted at in 
Remark 5.2. Also there is the surprising stochastic equation in (|5.8j) . There is 
also a related result given in Proposition 2.2.. One may wonder if properties of 
this sort only hold for beta random variables. We show in the next result, which 
was first obtained in that there is a considerable generalization which leads 
to some interesting stochastic equations. 

Proposition 6.1. Let 1Z = M T / a (Fn) and Q = M T / a (Fc)) denote Dirichlet 
means with parameters {t / a, R) and (t/oi, Q) where R is a non-negative random 
variable, and Q is a random variable taking values in [0, 1]. Equivalently, 

ft=M(*,i)ft + (l-% J i))ii; onde = %,i ) Q + (l-%,i ) )Q. (6.1) 

Which implies that Q = M T i a {Fo) takes it values in [0,1]. If Q is a constant 
then M T i q (-Fq) = Q. Then the following results hold. 

(i) n l ' a X a . T = M T (F XaR i /a ), that is, 

■R}l a X a , T = pr^^X^r + (1 - p Ttl )RV a X a . 
(ii) P a ,r{Q) is a Dirichlet mean with parameters (r, P a ,o(Q))> o,nd satisfies, 

p*,AQ) = 0t,iKA&) + C 1 - MPaAQ) ( 6 - 2 ) 

where P' aT {Q!) is equivalent only in distribution to P a ,r{Q)> but Q is the 
same variable. 

Proof. The first result follows by noting 

rE[log(l + XXaR 1 ^)} = (r/a)E[log(l + X a R)}, 

which gives the negative logarithm of C T (A; 7?. 1 / a A QjT ) = C T / a (X a ; TV). For the 
second, evaluate C T {\P a , T {Q)) conditional on Q and then notice similar to 
Proposition 5.2, that the transform of order r coincides with the exponential of 

logC T/Q ((l + A) Q - 1; Q) = -^Efofrg], 
It remains then to apply (|5.1|) . to get the second equality in (|6.3p . □ 
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Notice that setting R = £ CT and Q = p^ a we recover Propositions 2.2 and 
Proposition 5.2. Setting R = X$ for < 5 < 1, leads to the identity )T = 

X^", since it follows from known properties of stable random variables that 

Xy a X a = X a $. Furthermore, if one chooses Q := Q[u) such for each fixed u it 
satisfies (|6.1|) . and for < u < 1 it is an exchangeable bridge, that is a random 
cumulative distribution function, then P a ,r(Q(w)) identifies a coagulation oper- 
ation as described in Pitman (|8(j), Lemma 5.18), see also Bertoin Q). In partic- 
ular, one recovers Pitman's (|78l ) coagulation as follows. Setting Q = -Pg jT / Q (w), 

means that Q = P§fi(u), leading easily to, P a fi(Ps,o(u)) = P a s,o(u), which im- 
plies 

Pa,T(Ps,T/a(u)) = P a 6,r(u). 

We shall discuss other applications of Proposition 6.1 and related identities 
elsewhere. 

7. Subordinators and symmetric generalized Linnik laws and 
processes 

Related to the previous section we first briefly discuss some properties of their 
corresponding GGC subordinators which shows that they can be used in appli- 
cations where stable and exponentially tilted stable subordinators are used but 
are generally much more easy to handle from a practical point of view. This is 
especially true for the tilted stable subordinators, whose practical usage often 
requires non-trivial simulation techniques. 
Similar to the case of Xa,r, it follows that 

T a (t) = lr -R}l a X^ r = Xa,r^ 1/Q = J T M T (F XaR ya), 

and 

T Q (r) = 7t P q , t (Q) 

are GGC (rjIaR 1 / ) and GGC(t, P a fi(Q)) variables, respectively. Where we 
are suppressing the fact that both 1Z and Q depend on (a,r). In fact T q (t) 
and T q (t) are GGC subordinators varying in r > 0. Let S a (t) denote a positive 

stable subordinator such that S a (l) = S a , and let S a (t) denote the subordinator 
with 

- logE[e _A ^ (t) ] = i[(l + X) a - 1], 

so that S a (l) = S a is a random variable with density e _t f a (t). We note that 
S a is much more challenging to simulate than S a . 
It follows that, 

S a (-/z.K) = T a (t) and 5 a ( 7 xQ) £ f Q (r). 

It is evident that all such T a (r) involve X a ^ T , and one can work with Xa.r instead 
of S a . For T a (r) the extra randomization allows one to avoid sampling variables 
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such as S a . For instance, assuming one can sample <2, then first conditioning on 
Q = p one can sample the variables (7 T , P QjT (p), Q) However, in many cases it 
is not necessary to sample P a ^ T (p). For example, using Proposition 5.3, for any 
< cr < 1 



In this case one can work with 71 and a variable with density £l a ,<j- This also 
implies that, due to independent increments, one has an explicit description 
of the finite dimensional distributions. Of course if one can sample 1Z and Q 
directly these variables need not correspond to Dirichlet means. In this sense 
X a ^ T and P a , T (p) are basic to T a and T a . Naturally if 1Z and Q are Dirichlet 
means, this often leads to additional simplicities, and in many cases it is not 
necessary to sample Q or 1Z This is particularly true i f Q and R are not too 
complicated and < r < 1. One can also use James (|45l ) to obtain explicit 

densities. Lastly if lim Q ^o Q = Qo exists then by continuity properties of P a T , 



(((82), section 5.2), for fixed r, 

lim Q £ Qo, and lim f Q (r) = <y T Po, T (Qo)- 

Hence these Levy processes are attractive in terms of potential applications 
arising for instance in finance, financial econometrics or Bayesian statistics. With 
applications to finance in mind, it is quite natural to use these processes as 
Brownian time changes creating process B(T a (-)) and B(T a (-j), for B(-) an in- 
dependent Brownian motion, we take to have log characteristic function —A 2 . 
The characteristic functions of B(T a (r)) and B(r a (r)), can be expressed as 

Cz.{\ 2a -K) = e -^ A2Q ) and C ^{{\ + A 2 )" - 1; Q) = e -5*»« 1 +* , >'"- 1 >. 

Recall that, B(S a (-)) is a symmetric stable process of index (0, 2] and B(S a (-)), is 
a process that includes the NIG process when a = 1/2. When a = and Q = p, 
B(Tq(-)), is a variance-gamma (VG) process. The case of B(x a .g), corresponds 
to generalized Linnik processes considered by Pakes(j73j), see also (27), (61). We 
will focus on this case. 

It suffices to examine the random variables, 



where N is a standard Normal random variable. For general a and 9 > the 
extra randomization by N does not add much beyond our results for Xafi- 
However, when a < 1/2 we are able to obtain some interesting results which 



we describe below. In this case, we will use a result of Devroye (|26f ) which yields 
a tractable mixture representation for symmetric stable random variables of 
index between and 1. Devroye's (j26h result is not well known but as we shall 
show can be used to obtain a nice description of the density of B a (T a (6)) for all 
fixed 9 > 0. We do however stress that there are many applications requiring 
a > 1/2. 
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7.1. a < 12, Results based on Fejer-de la Vallee Poussin mixtures 



For symmetric stable random variables Ny/2S a for < a < 1/2, Devroye (|2 

shows that 

Ny/2S^=Y/Z^ 

where Y has a Fejer-de la Vallee Poussin density, 

2 



u){x) 



1 /sin(x/2) 



2tt\ x/2 



-oo < x < oo 



and Z = 7i (1 — ^2a) + 72^20 • It follows that the density of a symmetric stable 
of index between [0, 1] is, 



u(xy)y 2a e- y2a [(1 - 2a) + 2ay 2a ]dy. 



2a 



Hence as a mild extension of Devroye f (|26T ). Example B), that is the simple 
symmetric Linnik variable corresponding to 9 — a, we have for all 9 > 0, 



N^2x~ e = Y/W^ where W = ^(1 - £ 2a ) + ^ 2a 

7«. 7«. 

is a mixture of Pareto variables, having density for < a < 1/2 
0[(l + 20)w + (l-2cQ] 
a (l+ w ) 2+9 / Q ' 

Naturally this representation extends to all B(T a (-)), provided that a < 1/2. In 
particular, 

B(T a (6)) = Y/W^ where now W = -^-(1 - £ 2ct ) + -^T^a- 



72 



79.72. 



Quite interestingly the density of only requires information about the Laplace 
transform of 70 / a lZ. Let ip^ (x) and -0^ "* (x) denote the first and second deriva- 
tives of ipji(x). Then the density of W is given by 

= ^e-5**W[^g)(a;)(l - 2a) + x2a[(^\x)f^ - i/>%\x)]] 
From this, we close with an interesting identity, 

Proposition 7.1. For < a < 1/2, and 9 > 0, let V = /3i/2,i/2: / or 
—00 < a; < 00 



E 



2 -E 

7T 



4 Xo,» 



\y/2Xa,6 

Lj{xy)29y 2a [{\ + 29)y 2a + (1 - 2a)] 



(l + V 2a ) 



2a\2+e/a 
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Which is just the density of N y/lxafi- Additionally, for all fixed 9 > 0, the 
density of B(T a (9)) = N yJ2x a fi7l^ , satisfies 



E 



1 



2a 



u(xy)y' a n a fi(y 2a )dy 



Proof. The result follows from the derivation of the density described above 
using uj, in combination with a derivation of the density based on iV 2 = 271 V, 
and additionally Pitman and Yor f |85l), eq. (29)). □ 



8. General remark about rational case 



When a is rational it is known, see (|21l ; 1 9 It l96f). that S a is equivalent in distri- 



bution to a product of independent beta and gamma variables. Extending an 
argument in Chaumont and Yor ((21), p. 143-144) using the gamma duplication 
formula, it follows that for a = m/n for integers, m, n, such that m < n, and 
all 8 > —m/n, 




'm-X 

where all random variables are independent. Additionally 

)m tm— 1 

An implication of these relationships is that one may use the result of Springer 




and Thompson (88) to express their densities in terms of Meijer-G functions. 
In many cases these are equivalent to expressions in terms of generalized Gauss 
hypergeometric functions. Furthermore, it is known that Laplace transforms of 
Meijer-G functions are also Meijer-G functions, with known arguments. Hence, 
Proposition 3.1 and Proposition 5.1 show that in the rational case of a — 
m/n, one may express the generalized Mittag-LefHcr functions and densities 
for P a .e(p) in terms of Miejer-G functions. Such representations are not entirely 
appealing in many respects, for instance the density A m / n a is a much more de- 
sirable expression than its Miejer-G counterpart. However, from a computational 
viewpoint they are significant. This is due to the fact that Meijer-G functions, 
which constitute many special functions, are available as built-in functions in 
mathematical computational packages such as Mathematica or Maple. Natu- 
rally many of the quantities we discussed for general a can be expressed as the 
more general Fox-H functions However, in general, com- 



putations for these expressions are not yet available. Hence another contribution 
of our work is to give new explicit identities for a class of Meijer-G and Fox-H 
functions. That is to say quantities such as A a ^ give an explicit form to their 
corresponding Fox-H representation. We omit details of this representation, but 
it is not difficult to obtain. 
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9. An example from time series 

This last section is based on a question posed by Richard Davis which illustrates 
another case where additional randomization is helpful, albeit not obviously so. 



Davis and Resnick (|2a )[see also f (|18l ) .Theorem 13.3.1, p. 538-539.)], shows that 
in the infinite variance case, the sample correlation function of a moving average 
process converges to a variable of the form, 

B(S a ) £ N \fSg 

for all < a < 1, and where we omit various constants. The appearance of the 
square root seems to complicate matters. Looking at the square we get 

jl S'a d n . S a d 1 S a 



02 n / 2 C2 O 



where the last equality follows by noting that 471/2 = I/S1/21 & n d S^S 1/2 = 

S a /2- There are essentially two options left. Writing S a / 2 = = sH^Sa, 

one gets 

gq d_ A Q—l/a y A Al/a l/ot v 

~q — c C2 — D l/2 Aa ~ 4 '1/2 

Notice that the last two expressions consists of variables with known densities. 
So this illustrates another situation where how one randomizes matters. 



References 

[1] Aldous, D. and Pitman, J. (2006). Two recursive decompositions of 
Brownian bridge related to the asymptotics of random mappings. In memo- 
riam Paul-Andre Meyer: Seminaire de Probabilites XXXIX, 269-303, Lec- 
ture Notes in Math., 1874, Springer, Berlin. 

[2] Anh, V. V., Heyde, C. C. and Leonenko, N. N. (2002). Dynamic mod- 
els of long-memory processes driven by Levy noise. J. Appl. Probab. 39 
730-747. 

[3] Barlow, M., Pitman, J. and Yor, M. (1989). Une extension multidimen- 
sionnelle de la loi de l'arc sinus. In Seminaire de Probabilites XXIII (Azema, 
J., Meyer, P.-A. and Yor, M., Eds.), 294-314, Lecture Notes in Mathematics 
1372. Springer, Berlin. 

[4] Barndorff-Nielsen, O. E., Kent, J. and Sorensen, M.(1982) Normal 
variance-mean mixtures and z distributions. Internat. Statist. Rev. 50 145- 
159. 

[5] Barndorff-Nielsen, O. E., and Leonenko, N. N. (2005). Spectral 
properties of superpositions of Ornstein-Uhlenbeck type processes. Methodol. 
Comput. Appl. Probab. 7 (2005) 335-352. 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. James/Lamperti Laws 



34 



[6] Berestycki, J., Berestycki, N., and Schweinsberg, J. (2008) Small- 
time behavior of beta coalescents. Ann. Inst. H. Poincare Probab. Statist. 
44 214-238. 

[7] Bertoin, J. (1996). On the first exit time of a completely asymmetric stable 
process from a finite interval. Bull. London Math. Soc. 28 514-520. 

[8] Bertoin, J. (2006). Random fragmentation and coagulation processes, Cam- 
bridge University Press. 

[9] Bertoin, J. amd Le Gall, J.-F. (2006). Stochastic flows associated to 
coalescent processes. III. Limit theorems. Illinois J. Math. 50 147-181 

[10] Bertoin, J. and Goldschmidt, C. (2004). Dual random fragmentation 
and coagulation and an application to the genealogy of Yule processes. In 
Mathematics and computer science III: Algorithms, Trees, Combinatorics 
and Probabilities, M. Drmota, P. Flajolct, D. Gardy, B. Gittenberger (edi- 
tors), pp. 295-308. Trends Math., Birkhauser, Basel. 

[11] Bertoin, J. and Yor, M. (2001). On subordinators, self-similar Markov 
processes and some factorizations of the exponential variable. Electron. 
Comm. Probab. 6 95 106. 

[12] Bertoin, J. and Yor, M. (1996). Some independence results related to 
the arc-sine law. J. Theoret. Probab. 9 447-458. 

[13] Bertoin, J., Fujita, T., Roynette, B. and Yor, M. (2006). On a 
particular class of self-decomposable random variables : the duration of a 
Bessel excursion straddling an independent exponential time. Prob. Math. 
Stat, 26, p. 315-366. 

[14] Biane, P., Pitman, J., and Yor, M. (2001). Probability laws related 
to the Jacobi theta and Riemann zeta functions, and Brownian excursions. 
Bull. Amer. Math. Soc. (N.S.) 38 435-465. 

[15] Biane, P., and Yor, M. (1987). Valeurs principales associees aux temps 
locaux Brownicns. Bull. Sci. Math. (2) 111, 23 101. 

[16] BONDESSON, L. (1992). Generalized gamma convolutions and related 
classes of distributions and densities. Lecture Notes in Statistics, 76. 
Springer- Verlag, New York. 

[17] Bourgade, P., Fujita, T. and Yor, M. Euler's formulae for ((2n) and 
products of Cauchy variables. Electron. Comm. Probab. 12 73-80. 

[18] Brockwell, P. J., Davis, R. A. (2002) Introduction to time series and 
forecasting. Second edition. With 1 CD-ROM (Windows). Springer Texts in 
Statistics. Springer- Verlag, New York. 

[19] Carlton, M. A. (2002). A family of densities derived from the three- 
parameter Dirichlet process. J. Appl. Probab. 39, 764-774. 

[20] Chamati, H. and Tonchev, N. S. (2006). Generalized Mittag-LefHer 
functions in the theory of finite-size scaling for systems with strong 
anisotropy and/or long-range interaction. J. Phys. A 39 469-478. 

[21] CHAUMONT, L. AND YOR, M. (2003). Exercises in probability. A guided 
tour from measure theory to random processes, via conditioning. Cambridge 
Series in Statistical and Probabilistic Mathematics, 13, Cambridge Univer- 
sity Press. 

[22] Cifarelli, D. M. AND Regazzini, E.. (1979). Considerazioni generali 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. J ames / Lamperti Laws 



35 



sull'impostazione bayesiana di problemi non parametrici. Le medie asso- 
ciative nel contesto del processo aleatorio di Dirichlet I, II. Riv. Mat. Sci. 
Econom. Social 2 39-52 (1979). 

[23] Cifarelli, D.M. and Regazzini, E. (1990). Distribution functions of 
means of a Dirichlet process. Ann. Statist. 18, 429-442 (Correction in Ann. 
Statist. (1994) 22, 1633-1634). 

[24] Cifarelli, D.M. and Mellili, E. (2000). Some new results for Dirichlet 
priors. Ann. Statist. 28, 1390-1413. 

[25] DAVIS, R., Resnick, S. (1986). Limit theory for the sample covariance 
and correlation functions of moving averages. Ann. Statist. 14 533-558. 

[26] Devroye, L. (1984). Methods for generating random variates with Polya 
characteristic functions. Statist. Probab. Lett. 2 257-261. 

[27] Devroye, L. (1990). A note on Linnik's distribution Statist. Probab. Lett. 
9 305-306. 

[28] Devroye, L. (1996). Random variate generation in one line of code, in: 
1996 Winter Simulation Conference Proceedings, ed. J.M. Charnes, D.J. 
Morrice, D.T. Brunncr and J.J. Swain, pp. 265-272, ACM. 

[29] Dong, R., Martin, J. and GOLDSCHMIDT, C. (2005). Coagulation- 
fragmentation duality, Poisson-Dirichlet distributions and random recursive 
trees. Ann. Appl. Probab. 16 1733-1750. 

[30] DjRBASHlAN, M. M.(1993) Harmonic analysis and boundary value prob- 
lems in the complex domain. Translated from the manuscript by H. M. Jer- 
bashian and A. M. Jerbashian Operator Theory: Advances and Applications, 
65. Birkhauser Verlag, Basel. 

[31] Enriquez, N, Lucas, C, and Simenhaus, F. (2009). The Arcsine law 
as the limit of the internal DLA cluster generated by Sinai's walk. Preprint, 
larXiv:0903.4831k rl. 

[32] Ethier, S. N. and Griffiths, R. C. (1993). The transition function of 

a Fleming- Viot process Ann. Probab. 21 1571-1590. 
[33] Evans, S.N, and Ralph, P. L. (2009) Dynamics of the time to the most 

recent common ancestor in a large branching population. To appear Ann. 

Appl. Probab. 

[34] FERGUSON, T.S. (1973). A Bayesian analysis of some nonparametric prob- 
lems. Ann. Statist. 1, 209-230. 

[35] Fox, C. (1961). The G and H functions as symmetrical Fourier kernels. 
Trans. Amer. Math. Soc. 98 395-429. 

[36] Fujita, T. and YOR,M. (2006). An interpretation of the results of the 
BFRY paper in terms of certain means of Dirichlet processes. Preprint. 

[37] Gnedin, A. and Pitman, J. (2006). Exchangeable Gibbs Partitions and 
Stirling Triangles. J. Math. Sci. 138 5674-5685. 

[38] Haas, B., Miermont, G., Pitman, J. and Winkel, M. (2008). Con- 
tinuum tree asymptotics of discrete fragmentations and applications to phy- 
logenetic models. Ann. Probab. 36 1790-1837. 

[39] Haas, B., Pitman, J. and Winkel, M. (2008). Spinal partitions and 
invariance under re-rooting of continuum random trees Ann. Probab. 37 
1381-1411. 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. J ames / Lamperti Laws 



36 



[40] Hilfer, R., and Seybold, H. J. (2006). Computation of the general- 
ized Mittag-Lcfflcr function and its inverse in the complex plane. Integral 
Transforms Spec. Funct. 17 637-652. 

[41] Hjort, N.L. and Ongaro, A. (2005). Exact inference for random Dirich- 
let means. Stat. Inference Stock. Process., 8 227-254. 

[42] Hu, J. (2009). Modelling subordinated stochastic processes with Student's 
t and generalized secant hyperbolic Increments: empirical study of specula- 
tive energy markets. Available at SSRN: [http: / / ssrn.com/ abstract =1317893 

[43] Huillet, T. (2000) On Linnik's continuous-time random walks. J. Phys. 
A 33 2631-2652. 

[44] ISHWARAN, H. and James. L.F.. (2001). Gibbs sampling methods for stick- 
breaking priors. J. Amer. Statist. Assoc. 96, 161-173. 

[45] James, L.F.(2009). Dirichlet mean identities and laws of a class of subor- 
dinators. To appear Bernoulli. arXiv:0708.0614 

[46] James, L.F. (2006). Gamma tilting calculus for GGC and Dirich- 
let means via applications to Linnik processes and occupation time 
laws for randomly skewed Bessel processes and bridges. Available at 
|http://arxiv.org/abs/math.PR/0610218| 

[47] James, L. F., Lijoi, A. and Prunster, I. (2008). Distributions of func- 
tional of the two parameter Poisson-Dirichlet process. Ann. Appl. Probab. 
18, 521-551. 

[48] James, L.F., Roynette, B. and Yor, M. (2008). Generalized Gamma 
Convolutions, Dirichlet means, Thorin measures, with explicit examples. 
Probab. Surv. 5 346-415. 

[49] James, L.F. and Yor, M.(2006). Tilted stable subordinators, 
gamma time changes and occuation time of rays by Bessel spiders. 
|http://arxiv.org/abs/math.PR/0701049| 

[50] Jayakumar, K., AND Pillai, R. N. (1996). Characterization of Mittag- 
Leffier distribution. J. Appl. Statist. Sci. 4 77-82. 

[51] Kasahara, Y. (1975). Spectral theory of generalized second order differ- 
ential operators and its applications to Markov processes, Japan J. Math., 
1 67-84. 

[52] Kasahara, Y. and Yano, Y. (2005). On a generalized arc-sine law for 
one-dimensional diffusion processes, Osaka J. Math., 42 1-10. 

[53] Kawazu, K., and Watanabe, S. (1971). Branching processes with im- 
migration and related limit theorems. Teor. Verojatnost. i Primenen. 16 
34-51. 

[54] Keilson, J. and Wellner, J. A. (1978). Oscillating Brownian motion. 
J. Appl. Probability 15 300-310. 

[55] Kilbas, A. A., Saigo, M., Saxena, R. K. (2004). Generalized Mittag- 
Lcffler function and generalized fractional calculus operators. Integral Trans- 
forms Spec. Fund. 15 31-49. 

[56] Kilbas, A. A. and Saigo, M. (2004). H-Transforms. Theory and Appli- 
cations, Chapman and Hall/CRC, Boca Raton, FL. 

[57] Kiryakova, V. (1994). Generalized fractional calculus and applications. 
Pitman Research Notes in Mathematics Series, 301. Longman Scientific and 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. James/Lamperti Laws 



37 



Technical, Harlow; copublished in the United States with John Wiley and 

Sons, Inc., New York. 
[58] Kiryakova, V. (1997). All the special functions are fractional differintc- 

grals of elementary functions. J. Phys. A 30 14, 5085-5103. 
[59] Kotani, S. and Watanabe, S. (1982). Krein's spectral theory of strings 

and generalized diffusion processes. Functional Analysis in Markov processes, 

235-259, Lecture Notes in Math., 923, Springer, Berlin- Heidelberg-New 

York. 

[60] Kotlarski, I. (1965). On pairs of independent random variables whose 
product follows the gamma distribution. Biometrika 52 289-294. 

[61] KozuBOWSKi, T. J. (2001). Fractional moment estimation of Linnik and 
Mittag-Lefner parameters. Stable non-Gaussian models in finance and econo- 
metrics. Math. Corn-put. Modelling 34, no. 9-11, 1023-1035. 

[62] Kyprianou, A. E., Rivero, V. (2008). Special, conjugate and complete 
scale functions for spectrally negative Lvy processes. Electron. J. Probab. 13 
57 1672-1701. 

[63] Kyprianou, A. E., and Pardo, J. C. (2008). Continuous-state branch- 
ing processes and self-similarity. J. Appl. Probab. 45 1140-1160. 

[64] Lageras, Andreas Nordvall; Martin-Lof, A. (2006) Genealogy for 
supercritical branching processes. J. Appl. Probab. 43 1066-1076. 

[65] Lamperti, J. (1967). Limiting distributions for branching processes. 
1967Proc. Fifth Berkeley Sympos. Math. Statist, and Probability (Berke- 
ley, Calif., 1965/66), Vol. II: Contributions to Probability Theory, Part 2 
pp. 225-241 Univ. California Press, Berkeley, Calif. 

[66] Lamperti, J. (1958). An occupation time theorem for a class of stochastic 
processes. Trans. Amer. Math. Soc. 88, 380-387. 

[67] Levy, P. (1939). Sur certains processus stochastiques homogenes. Compo- 
sitio Math. 7, 283-339. 

[68] Lin, G. D. (2001). A note on the characterization of positive Linnik laws. 
Aust. N. Z. J. Stat. 43 1, 17-20. 

[69] Lukacs, E. (1955) A characterization of the gamma distribution. Ann. 
Math. Statist. 26 319-324. 

[70] Mainardi, F., Pagnini, G., and Saxena, R. K. (2005). FoxH functions 
in fractional diffusion. J. Comput. Appl. Math. 178 321-331. 

[71] Mainardi, F., Luchko, Y., and Pagnini, G. (2001). The fundamental 
solution of the space-time fractional diffusion equation. Fract. Calc. Appl. 
Anal. 4 153-192. 

[72] Nagaev, S. V., and Wachtel, V. (2007) The critical Galton- Watson 

process without further power moments. J. Appl. Probab. 44 753-769. 
[73] Pakes, A. G. (1998) Mixture representations for symmetric generalized 

Linnik laws. (English summary) Statist. Probab. Lett. 37 213-221. 
[74] Patie, P. (2009). Exponential functional of a new family of Levy processes 

and self-similar continuous state branching processes with immigration, Bull. 

Sci. Math., 133 355-382. 
[75] Perman, M., Pitman, J. and Yor, M. (1992). Size-biased sampling of 

Poisson point processes and excursions. Probab. Theory Related Fields. 92, 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. James/Lamperti Laws 



: J ,8 



21-39. 

[76] Pillai, R. N. (1990). Mittag-Lcfflcr functions and related distributions. 

Ann. Inst. Statist. Math 42 157161. 
[77] Pitman, J. (1996). Some developments of the Blackwell-MacQueen urn 

scheme. Statistics, probability and game theory, 245-267, IMS Lecture Notes 

Monogr. Ser., 30, Inst. Math. Statist., Hayward, CA. 
[78] Pitman, J. (1999). Coalescents with multiple collisions. Ann. Probab. 27 

1870-1902. 

[79] Pitman, J. (2003). Poisson-Kingman partitions. In Science and Statis- 
tics: A Festschrift for Terry Speed. (D.R. Goldstein, Ed.), 1-34, Institute of 
Mathematical Statistics Hayward, California. 

[80] Pitman, J. (2006). Combinatorial stochastic processes. Lectures from the 
32nd Summer School on Probability Theory held in Saint-Flour, July 7-24, 
2002. With a foreword by Jean Picard. Lecture Notes in Mathematics, 1875. 
Springer- Verlag, Berlin. 

[81] Pitman, J., and Yor, M.(1992). Arcsine laws and interval partitions 
derived from a stable subordinator. Proc. London Math. Soc. 65 326-356. 

[82] Pitman, J. and Yor, M.. (1997) The two-parameter Poisson-Dirichlet 
distribution derived from a stable subordinator. Ann. Probab. 25, 855-900. 

[83] Pitman, J., and Yor, M. (1997) On the relative lengths of excursions 
derived from a stable subordinator. In: Seminaire de Probabilites XXXI 
(Azema, J., Emery, M. and Yor, M., Eds.), 287-305, Lecture Notes in Math- 
ematics 1655. Springer, Berlin. 

[84] Pitman, J. and Yor, M.. (2003). Infinitely divisible laws associated with 
hyperbolic functions, Canad. J. Math. 55 292-330. 

[85] Pitman, J., and Yor, M.(2004). Some properties of the arc-sine law 
related to its invariance under a family of rational maps. A festschrift for 
Herman Rubin, 126-137, IMS Lecture Notes Monogr. Ser., 45, Inst. Math. 
Statist., Beachwood, OH. 

[86] Prabhakar, T. R. (1970). On a set of polynomials suggested by Laguerre 
polynomials. Pacific J. Math. 35 213-219. 

[87] Slack, R. S. (1968). A branching process with mean one and possibly 
infinite variance. Z. Wahrscheinlichkeitstheorie und Verw. Gebiete 9 139- 
145. 

[88] Springer, M. D., and Thompson, W. E.(1970). The distribution of 
products of beta, gamma and Gaussian random variables. SIAM J. Appl. 
Math. 18 721-737. 

[89] Teh, Y.W. (2006). A Hierarchical Bayesian Language Model based on 
Pitman- Yor Processes. Coling/ACL 2006. 

[90] Watanabe, S. (1993). Generalized arc-sine laws for one-dimensional dif- 
fusion processes and random walks, Stochastic analysis (Ithaca, NY, 1993), 
157-172, Proc. Sympos. Pure Math., 57, Amer. Math. Soc, Providence, RI, 
1995. 

[91] Williams, E. J. (1977) Some representations of stable random variables 

as products. Biometrika 64, 167-169. 
[92] Yaglom, A. M. (1947). Certain limit theorems of the theory of branching 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



Lancelot F. James/Lamperti Laws 



39 



random processes. (Russian) Doklady Akad. Nauk SSSR (N.S.) 56 795-798. 

[93] Yano, Y. (2006). On the occupation time on the half line of pinned diffu- 
sion processes. Publ. Res. Inst. Math. Set. 42 (2006), 787-802. 

[94] Yano, K. and Yano, Y. (2008). Remarks on the density of the law of the 
occupation time for Bessel bridges and stable excursions. Statist. Probab. 
Lett, 78 2175-2180. 

[95] K. Yano, Y. Yano and M. Yor. (2009). On the laws of first hitting times 
of points for one-dimensional symmetric stable Levy processes. Seminaire de 
Probability XLII, 187-227, Lecture Notes in Math., 1979, Springer, Berlin, 
2009. 

[96] Zolotarev, V. M. (1957) Mellin-Stieltjes transformations in probability 
theory. (Russian) Teor. Veroyatnost. i Primenen. 2 444-469. 

[97] Zolotarev, V. M. (1957) More exact statements of several theorems in 
the theory of branching processes. Teor. Veroyatnost. i Primenen. 2 256- 
266. 



imsart-generic ver. 2007/02/20 file: ArLamperti.tex date: August 31, 2009 



